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ABSTRACT 

Spectrum sensing is an essential functionality that enables 
cognitive radios to detect spectral holes and opportunistically 
use under-utilized frequency bands without causing harmful 
interference to primary networks. Since individual cognitive 
radios might not be able to reliably detect weak primary sig- 
nals due to channel fading/shadowing, this paper proposes a 
cooperative wideband spectrum sensing scheme, referred to 
as spatial-spectral joint detection, which is based on a linear 
combination of the local statistics from spatially distributed 
multiple cognitive radios. The cooperative sensing problem 
is formulated into an optimization problem, for which subop- 
timal but efficient solutions can be obtained through mathe- 
matical transformation under practical conditions. 

Index Terms — Spectrum sensing, distributed detection, 
nonlinear optimization, and cognitive radio. 

1. INTRODUCTION 

As an essential functionality of cognitive radio (CR) networks 
ITl, spectrum sensing needs to reliably detect weak primary 
radio signals of possibly-unknown formats. Generally, spec- 
trum sensing techniques can be classified into three categories: 
energy detection, matched filter coherent detection |2), and 
cyclostationary feature detection. Since non-coherent energy 
detection is simple and able to generate the spectrum-occupancy 
information quickly, we adopt it as the building block for con- 
structing the proposed wideband spectrum sensing schemes. 

The literature on wideband spectrum sensing for CR net- 
works is limited. An earlier approach is to use a tunable nar- 
rowband bandpass filter at the RF front-end to sense one nar- 
row frequency band at a time, over which the existing narrow- 
band spectrum sensing techniques can be applied. In order to 
search over multiple frequency bands at a time, the RF front- 
end needs a wideband architecture, and spectrum sensing usu- 
ally operates over an estimate of the power spectral density 
(PSD) of the wideband signal. In [31, wavelet transformation 
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was used to estimate the PSD over a wide frequency range 
given its multi-resolution features. However, no prior work 
attempts to make decisions over multiple frequency bands 
jointly, which is essential for implementing efficient CR net- 
works. 

In this paper, we consider the situation in which spectrum 
sensing is compromised by destructive channel conditions be- 
tween the target-under-detection and the detecting cognitive 
radios, which makes it hard to distinguish between a white 
spectrum and a weak signal. We propose a cooperative wide- 
band spectrum sensing scheme that exploits the spatial diver- 
sity among cognitive radios to improve the sensing reliabil- 
ity. The cooperation is based on a linear combination of local 
statistics from spatially distributed cognitive radios [4| [5!|, 
where the signals are assigned different weights according to 
their positive contributions to joint sensing. In such a sce- 
nario, we treat the design of distributed wideband spectrum 
sensing as a spatial-spectral joint detection problem, which 
is further formulated into an optimization problem with the 
objective of maximizing the overall opportunistic through- 
put under constraints on the interference to primary users. 
Through mathematical reformulation, we derive a suboptimal 
but efficient solution for the optimization problem, which can 
considerably improve sensing performance. 

2. SYSTEM MODEL 

Consider a primary communication system (e.g., multicarrier 
based) over a wideband channel that is divided into K non- 
overlapping subchannels. At a particular time, some of the 
K subchannels might not be used by the primary users and 
are available for opportunistic spectrum access. Multiuser or- 
thogonal frequency division multiplexing (OFDM) schemes 
are suitable candidates for such a scenario since they make it 
convenient to nullify or activate some portion of multiple nar- 
row bands. We model the detection problem over the subband 
k as one to choose between hypothesis Tio.fc ("0")> which rep- 
resents the absence of primary signals, and hypothesis Tii^k 
("1"), which represents the presence of primary signals. An 
illustration where only some of the K bands are occupied by 
primary users is illustrated in Fig. [T] The crucial task of spec- 
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Fig. 1. Illustration of the occupancy of a multiband channel. 

trum sensing is to sense the K frequency bands and identify 
spectral holes for opportunistic use. For simplicity, we as- 
sume that the high-layer protocols guarantee that all CRs keep 
quiet during the detection such that the main spectral power 
under detection is emitted by the primary users. 

Consider a multi-path fading environment, where h{l), 
Z = 0, 1, . . . , L — 1, denotes the discrete-time channel impulse 
response between the primary transmitter and a CR receiver 
with L equal to the number of resolvable paths. The received 
baseband signal at the CR front-end can be expressed as 

r{n) ^^h{l)s{n- l)+v{n), n = 0, 1, . . . , A^o-l (1) 

1=0 

where s{n) represents the primary transmitted signal with 
cyclic prefix at time n and v{n) is additive complex white 
Gaussian noise with zero mean and variance cr^, i.e., v{n) ^ 
CJV (O, CTy). In a multi-path fading environment, the wide- 
band channel exhibits frequency-selective features and its dis- 
crete frequency responses are given by 

L-l 

^'^ = ^ E Hn)e'^'-"''^''\ fc - 0, 1, . . . , A'-l (2) 

where L < Nq. We assume that the channel is slowly varying 
and the channel frequency responses {-fffcj^Jo^ ^^^y 
much during a detection interval. In the frequency domain, 
the received signal at each subchannel can be estimated by 
computing its discrete Fourier transform (DFT): 



No-l 



n=0 



where Sk is the primary signal at subchannel k and 



Vk = 



L-l 

E 

n=0 



v{n)e 



-j27rnk/No 



(3) 



(4) 



is the received noise in the frequency domain. Note that Vk ~ 
CJV (0, CT^) since v{n) - CM (O, cr^) and the DFT is a uni- 
tary linear operation. Without loss of generality, we assume 
that the transmitted signal Sk, channel gain Hk, and additive 
noise Vk are independent of each other 

To decide whether the fc-th subchannel is occupied or not, 
we test the following binary hypotheses: 

Ho,fc: Rk = Vk, A: = 0,1,..., if -1 

ni,k : Rk = HkSk + Vfe, A: = 0, 1, . . . , X - 1 (5) 

For each subchannel fc, we compute the test statistic as the 
sum of received energy over an interval of M samples, i.e.. 



Yk = 



M-l 

E i^fcMi' 

m=0 



0,1, 



.K-l 



(6) 



and the decision rule is given by 

'Hi,k 

Yk I Ik (7) 

where 7^ is the corresponding decision threshold. For sim- 
plicity, we assume that the transmitted signal at each subchan- 
nel has unit power, i.e., E (|S'/£p) = 1; this assumption holds 
when the primary radios adopt uniform power transmission 
strategies given no channel knowledge at the transmitter side. 

According to the central limit theorem for large M, Yk is 
asymptotically normally distributed with mean 



iE(n) 

and variance 

Var(rfc) = 



Mai 
M{al 

2Ma* 
2M (al - 



\Hk 



2\Hk 



Ti-ik 



Ti-o.k 
'Hi,k 



(8) 



(9) 



for k ~ 0,1, . . . , K ~ 1. Thus, assuming large M, we have 
Yk^Af{E{Yk),Yi,r{Yk)). 

Using the decision rule in (|7|, the probabilities of false 
alarm and detection at subchannel k can be respectively cal- 
culated approximately as 



^fH7fc) =Pr(i"fc >7fc|Ho) 



7fc 



Mat 



(10) 



and 



(11) 



/ j,-M{al + \Hk\') \ 

\a,^2M{al + 2\Hkn) 
where Q denotes the tail probability of the standard normal 
distribution. The choice of threshold jk leads to a tradeoff 
between the probabilities of false alarm and miss Pm = 1 — 
Pd- Specifically, a higher threshold will result in a smaller 
probability of false alarm but a larger probability of miss, and 
vice versa. 

3. SPATIAL-SPECTRAL JOINT DETECTION 

Suppose that N spatially distributed cognitive radios collabo- 
ratively sense a wide frequency band. By combining the local 
statistics from individual cognitive radios at the fusion cen- 
ter, which can be one of the CRs, the network can make a 
better decision on the presence or absence of primary signals 
on each of the K subchannels. The cooperation assumes a 
separate control channel, through which the statistics of in- 
dividual CRs are transmitted to the fusion center Let Yfe(n) 
denote the received energy in the fc-th subchannel at cognitive 
radio n. For each subchannel, these statistics can be written 
in a vector as Yk = [Yk{0), Yk{l), . . . ,YkiN - l)f. 

To exploit the spatial diversity, we linearly combine the 
summary statistics from spatially distributed cognitive radios 
at each subchannel fc to obtain a final test statistic: 



N-1 



Zk 



'Wk{n)Yk{n) = wlYk 



where KJfe = [■Wk{0),Wk{l), ■ ■ ■ ,Wk{N — I)] ai-e the com- 
bining coefficients for subchannel k, which can be compactly 
written in a matrix as W = [wo Wi . . . wk-i]- Note that 
Wk {n) > 0, for every k and n. 

Since the entries in Yk are normally distributed, the test 
statistics {zk}^^^ are also normally distributed with means 

where 1 is an all-one vector, and variances 



Var (zk) 



2MatwTwi 



n 



0,k 



[ 2Malwl [all + 2diag{Gk)]wk Hi^k 

(14) 

where Gfc = [|iJfc(0)|MiIfc(l)p, . . . , |iIfc(7V - 1)|2]^ are 
the squared magnitudes of the channel gains between the pri- 
mary transmitter and the N CR receivers for subchannel k. 

In order to decide the presence or absence of the primary 
signal in subchannel k, we use the following binary test 

Ti-i.k 

Zk I Ik, k = 0,l,...,K-l. (15) 
'Ho,k 

Accordingly, the detection performance in terms of the prob- 
abihties of false alarm and detection are given by 



Ik 



alJ2Mwlwk 



and 



P^''\wk,jk)^Q 



(16) 



Ik 



Mwl {all 



Gk) 



a^j2Mwl [all + 2diag(Gfe)] Wk 



(17) 

For compactness of notation, we collect the probabilities of 
false alarm and detection over the K subchannels into vectors 
P/(W,7) and Pd{W Thus, the probabilities of miss 
can be represented as Pm(W^,7) = 1 — Pd(W^,7). 

Our goal is to maximize the opportunistic rate while meet- 
ing some constraints on the interference to the primary com- 
munication system. Let Vk denote the throughput achiev- 
able over the fc-th subchannel if used by cognitive radios, and 
r = [ro, ri, . . . , rx-i]'' ■ Since 1 — pj*^^ measures the oppor- 
tunistic spectrum utilization of subchannel fc, we define the 
aggregate opportunistic throughput capacity as 



i?(W,7) = r^ [1-P/(W,7)] 



(18) 



For a widband primary system, the impact of interfer- 
ence induced by cognitive devices can be characterized by 
a relative priority vector over the K subchannels, i.e., c = 
[cq, ci, . . . , ck-i]^, where Ck indicates the cost incuiTed if 
the primary user at subchannel k is interfered with. As such, 
we define the aggregate interference to the primary user as 



Pm{W Consequently, the spatial-spectral joint de- 



tection problem is formulated as 



max R (W, 7) 



(PI) 



(19) 
(20) 
(21) 



S.t. C^PmiW,j) <S 

Pf{W,j)^f3 

where a = [ao, aR-iY" and (3 = [/3o, • ■ • , pk-\ 

Finding the exact solution for the above problem is diffi- 
cult since for any fc, P^'^\wk, jk) and p'^\wk-, Jk) are nei- 
ther convex nor concave functions according to ( fTSI l and ( fTTI l. 
To jointly optimize W and 7, we can show that (PI) can be 
reformulated into an equivalent form with convex constraints 
and an objective function lower bounded by a concave func- 
tion under the following conditions: 



< afc < i and < /3a; < ^, k = 0, . 



.,K-l. (22) 



Through maximizing the lower bound of the objective func- 
tion, we are able to obtain a good approximation to the opti- 
mal solution of the original problem. 

First, we show how to transform the nonconvex constraints 
in ( I20I ) and (l2Tl i into convex constraints by exploiting the 
monotonicity of the Q-function. Substituting ( fT6l l into the 
constraint ( |2TI ). we have 



Q-\l3k)\l2Mwlwk < 5 - Mwll 



(23) 



where Q^^{Pk) > given < 1/2. From ([T7), the con- 
straint ( |20b can be expressed as 

2M<[a?/ + 2diag(G.)]«,. < " ^^^'^ + ^'^^ 



o-vQ ^{l-ak) 

(24) 

given Qffc < 1/2 and Q ^(1 — a/c) < 0. Since the left-hand 
side on the constraint ( |23] l is convex and the right hand side 
is linear in (7^., Wk), (|23] | defines a convex set for (7^,, Wk). 
Similarly, ( |24] | is also a convex constraint. 

Then, we reformulate (PI) by introducing a new variable 



fik - (Jv^2Mwl [all + 2diag(Gfc)] Wk. (25) 

By defining 7^ = ^k/ ^'k and w'^. — Wk/nk, the constraints 
(|23] | and (l24l i can be further written as 



Q-\(3k)^J2Mw'^' < < ^ - Ml^w'f, (26) 



and 



j', - M {all + Gk) < <a,Q-i(l-«fe). (27) 
Note that ( |27] l is actually a linear constraint in (7^, lo'^,). The 
constraint ( fT9l ) now becomes 



l^c - CkQ [I'k ~ M {all + Gk) wij < e, (28) 

fe=0 

which can be shown to be convex by the following result. 



Lemma 1 If •yl < Af (c^l + Gfe)'^ wj,, then the function 

Q ( 7fc - M [all + w'A is concave in {7^, w',^}. 

By changing the variables W' and 7', P^- (10^,7^) can 
be expressed as 




2<Miag(Gfc)t^;, 



(29) 



From the the Rayleigh-Ritz theorem fSl, we have 



min|i/fc(n)| < 



w'k diag(Gfe)u;'fc 



< max I iJfc(n)r (30) 



Define a new function 



0-2 



2min|iJfc(n)|2 

n 

(31) 



which can be shown to be convex by the following result. 

Lemma 2 Ifjl > (tIAII^w'^., then the function gkil'k, w'/,) 
is convex in {7^, ii;'^,}. 

Since P^''^ {wk, jk) < dkilk^'^'k)^ '^he objective function 

in (PI) can be lower bounded by ^^Sq rk [1 — gk ilk^'^'k)]' 
which is a concave function. Thus, an efficient suboptimal 
method to solve (PI) is to maximize the lower bound of its 
objective function, i.e.. 



max 
w'n' 



St. 



AT-l 

E 

k=0 
K-1 

-E 

fc=0 



rk [1 -5fc [ikiw'k)] 



(P2) 



CfcC 



l^c 



Q-\(5k)^Iw' 



I T 



ik 



M {all + Gk) w'k < (JvQ'^l - ak). 

Implied by the practical conditions in ( l22l i. this problem is a 
convex optimization problem and can be solved efficiently. 

4. SIMULATIONS 

Suppose that two CRs cooperatively sense a multiband OFDM 
system with 8 subbands. For each subband, it is expected 
that the opportunistic spectrum utilization is at least 50%, 
i.e., Pk = 0.5, and the probability that the primary user is 
interfered is at most ak = 0.1. It is assumed that al — 1 
and M = 100. Other parameters are given in Table [T] Fig. 
12] shows result of solving (P2), which maximizes the oppor- 
tunistic throughput subject to the constraints on the interfer- 
ence. We observe that the joint detection results in much 
higher opportunistic throughput than the algorithms without 
cooperation. Note that the increase in the throughput of the 
joint optimization scheme becomes rather slow as we relax 
the interference constraint because the interaction between 7 
and W pushes the system to an operating point at which the 
throughput is more limited by /3 than by e. 



Table 1. Parameters used in simulations 
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Fig. 2. Aggregate opportunistic throughput capacity vs. the con- 
straint on the aggregate induced interference. 

5. CONCLUSION 

In this paper, we have proposed a spatial-spectral joint de- 
tection framework for distributed wideband spectrum sens- 
ing in cognitive radio networks, within which the cooperation 
among spatially distributed cognitive radios is optimized over 
multiple frequency bands. By exploiting the inherent struc- 
ture of the formulation, we have developed suboptimal but 
efficient solutions for the non-convex optimization problem. 
This paper establishes important principles for the design of 
distributed wideband spectrum sensing algorithms in cogni- 
tive radio networks. 
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